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We study the equilibrium correlations of a Bose gas in an elongated three-dimensional harmonic 
trap using a grand-canonical classical-field method. We focus in particular on the crossover from a 
quasicondensate regime characterized by strong phase fluctuations to a so-called true condensate that 
is phase coherent across its spatial extent. Choosing realistic experimental parameters, we quantify 
the density fluctuations and phase coherence of the atomic field as functions of the system temper- 
ature. We identify the onset of Bose condensation through analysis of both the generalized Binder 
cumulant appropriate to the inhomogeneous system, and the suppression of the effective many-body 
T matrix that characterizes interactions between condensate atoms in the finite-temperature field. 
We find that the system undergoes a second-order transition to condensation near the critical tem- 
perature for an ideal gas in the strongly anisotropic three-dimensional geometry, but remains in a 
strongly phase-fluctuating quasicondensate regime until significantly lower temperatures. We char- 
acterize the crossover from a quasicondensate to a true condensate by a qualitative change in the 
form of the non-local first-order coherence function of the field, and compare our results to those of 
previous works employing a density-phase Bogoliubov-de Gennes analysis. 
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I. INTRODUCTION 

Advances in the experimental control, observation and 
manipulation of quantum-degenerate dilute atomic gases 
have led to a large body of work focusing on the role 
of geometry and dimensionality in the physics of these 
quantum fluids [TT41 1J . In homogeneous systems, Bosc- 
Einstein condensation (BEC) of an ideal gas is precluded 
in dimensions d < 3, and long-range order in interacting 
systems is prohibited in such low dimensionalities |12j . 
Although some rigorous results exist for the harmonic 
trapping geometries typical of experimental dilute-gas 
systems [T3], these systems are theoretically less well 
characterized than their homogeneous counterparts. Ex- 
perimentally, these systems may exhibit quasicondensate 
behavior, characterized by large phase fluctuations and 
comparatively subdued density fluctuations, similar to 
that predicted to occur in homogeneous low-dimensional 
systems [TH [15] . On the other hand, the finite size of such 
systems can induce some phase coherence across the spa- 
tial extent of the atomic sample. This might loosely be 
associated with "finite-size" condensation, which will not 
exhibit the extensivity property of a formal, thermody- 
namic Bose condensate 16J. The physics of finite-sized 
and inhomogeneous quantum liquids in low dimensions 
are often more subtle than the infinite homogeneous sys- 
tems for which rigorous results are known (see, for ex- 
ample, Refs. [T7] HH]), and the effects of finite-size con- 
densation and its relationship to quasicondensation and 
superfluidity in low-dimensional systems remain in gen- 
eral somewhat unclear. 



A low-dimensional system of particular interest is 
the weakly interacting Bose gas in an elongated (cigar- 
shaped) three-dimensional (3D) harmonic trap [TJ ® 0]- 
El M ED [MM]- In the limit of an extreme trap 
anisotropy, such that the oscillator energy spacing in 
the two tightly confined (transverse) dimensions is much 
larger than the energy scales associated with interactions 
and thermal fluctuations, this system can be regarded 
as truly one-dimensional [2SH2Z1- Using a Bogoliubov- 
de Gennes (BdG) approach in density and phase fluctu- 
ations, Petrov et al. [35] determined the phase diagram 
for this system, and identified a regime of quasiconden- 
sate behavior. Their analysis nevertheless showed that, 
at sufficiently low temperatures, long- wavelength fluctu- 
ations of the phase are suppressed. The system is then 
phase coherent across a large spatial extent, and is said 
to contain a "true" condensate. 

Importantly, Bose gases in less severely elongated har- 
monic traps, in which the energy scales associated with 
interactions and thermal fluctuations are not smaller 
than the transverse oscillator spacings (and which are 
therefore formally 3D), can also exhibit quasiconden- 
sate behavior. This was demonstrated theoretically by 
Petrov et al. [25], who adapted their BdG density-phase 
approach to the elongated-3D geometry, and has been 
observed experimentally by a number of groups [191 1211 - 
I24j . As in the purely one-dimensional case, phase fluctu- 
ations in this system become suppressed at low temper- 
atures, yielding a gradual crossover from a quasiconden- 
sate regime to a true condensate. The pronounced effect 
that interactions have on the behavior of this system is 
seemingly in stark contrast to the familiar case of unam- 
biguously three-dimensional harmonically trapped Bose 
gases, in which interactions serve only to slightly decrease 
the critical temperature for BEC (see, e.g., Ref. [29"]). 
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As theoretical techniques based on a decomposition 
into density and phase fluctuations EE] ESI [30] make 
no reference to an underlying condensate, the nature of 
the transition to condensation is obscured in such ap- 
proaches. In particular, it is not clear from such calcu- 
lations to what extent the familiar picture of BEC as a 
second-order transition remains relevant to the physics 
of the elongated case. A thorough understanding of the 
relationship between condensate and quasicondensate in 
this straightforward scenario would seem to be a natural 
starting point for understanding the role of (potentially 
strictly finite-size) condensation in low-dimensional sys- 
tems. 

In this article we investigate the elongated 3D harmon- 
ically confined weakly interacting Bose gas using a grand- 
canonical classical-field method 31 33 that allows us to 
carefully characterize the condensate and quasiconden- 
sate, and the relationship between the two. Although 
related methods have been used in several studies of 
the condensate-quasicondensate crossover in (quasi-)one- 
dimensional systems [3~4H3"5] , the only prior classical-field 
investigation of the elongated 3D system is that of Kadio 
et al. [40] . The authors of Ref . [40] calculated the phase- 
coherence length of the elongated 3D system, and found 
approximate agreement with the predictions of Petrov et 
al. EE]. However, their study employed somewhat ad 
hoc techniques — based on ideal-gas arguments — to es- 
timate the temperature of the strongly fluctuating equi- 
librium state of the field, and was not able to access the 
statistics of the condensate mode itself, due to technical 
limitations [40] . As such, no detailed, quantitative ex- 
position of the relation between condensation and quasi- 
condensation in the elongated 3D system is available in 
the previously published literature. 

Here we present a detailed, quantitative analysis of the 
quasicondensate regime of the elongated 3D gas, focusing 
in particular on the nature of condensation in the system. 
We characterize condensation in the system by examining 
the number fluctuations of the orbital associated with the 
largest eigenvalue of the one-body density matrix, and 
the anomalous correlations of the part of the field orthog- 
onal to this orbital, and show that the gas exhibits mean- 
ingful Bose condensation in the quasicondensate regime. 
Moreover, we find that the onset of quasicondensation in 
the system is accompanied by a minimum in the effec- 
tive interaction strength between condensate atoms that 
is consistent with the system undergoing a second-order 
transition to condensation. We find that this occurs at 
a temperature slightly below the transition temperature 
of the corresponding ideal-gas model, and that the ap- 
pearance of a quasicondensate in the system is there- 
fore associated with the condensation transition. The 
crossover from a quasicondensate to a true condensate 
at lower temperatures can thus be understood in terms 
of the correlations of the complementary noncondensed 
component of the field, which we find are inconsistent 
with well-defined quasi-particles [i.e., with Gaussian, or 
Hartree-Fock-Bogoliubov (HFB) correlations [HJ [35]] & t 



high temperatures, but come to be more consistent with 
HFB correlations as the system temperature is reduced. 

This article is organized as follows: In Sec. [H] we 
describe the theoretical methods used in our analysis. 
We briefly explain our classical-field approach and the 
stochastic projected Gross-Pitaevskii equation (SPGPE) 
with which we describe the low-energy region of the sys- 
tem (Sec. II A), and discuss how we use it to calculate 



observables of interest (Sec. II B). In Sec. Ill we de- 



fine the physical parameters of the system we investi 



gate (Sec. Ill A I and present our analysis of its physical 



properties at varying temperature and fixed total atom 
number (Sees. IIIB - HIE I. In Sec. IV we summarize our 



results and present our conclusions. 

II. THEORETICAL METHODS 
A. Stochastic projected Gross-Pitaevskii equation 

The stochastic projected Gross-Pitaevskii equation 
(SPGPE) method, developed in Refs. [3TJ S3 HI], has 
been reviewed in detail, together with other projected 
classical-field methods, in Ref. [32] (see also Refs. [33l 
US]). For the reader's convenience, we briefly describe 
the relevant details of the formalism here. 

Formally, the physics of the harmonically trapped di- 
lute Bose gas is governed by the second-quantized Hamil- 
tonian 

H= fdr¥(r)H sp ${r) (1) 
+ \ j dr j dr ' ^(^VMr - r')*(r')*(r), 



where the single-particle Hamiltonian is 
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and U(r) is the exact interatomic potential. We intro- 
duce a single-particle subspace L spanned by eigenmodes 
Y n (r) of the single-particle Hamiltonian [H sp Y n (r) = 
e n Y n (r)} with energies e„ less than a cutoff energy -E max , 
and a complementary subspace comprised of the remain- 
ing high-energy modes. Provided £ max is chosen such 
that the high-energy modes are essentially unoccupied, 
the dynamics of these modes can be integrated out to ob- 
tain an effective Hamiltonian for the low-energy (coarse- 
grained) Bose field = SneL ^n^.( r ). as shown 
by Morgan [IS]. Atomic interactions described by the 
effective Hamiltonian are mediated by an approximate 
two-body T matrix, and the interaction can thus be rig- 
orously approximated by a "contact" potential, with a 
renormalized coupling constant Uq. In practice the cor- 
rection due to the finite momentum cutoff is small [47] . 
and so we assume the standard s-wave coupling constant 
Uq = 4irh 2 a/m, with a the s-wave scattering length. The 
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low-energy Hamiltonian then takes the form 



H h = /rfr*i(r)F sp * L (r) 



U 



dr^(r)*[(r)* L (r)* L (r), (3) 



which defines an effective field theory [48 for the coarse- 
grained field ^L(r). 

We then further divide the low-energy region L into 
a coherent region (or condensate band) C = {n : e n < 
£cut}, spanned by single-particle eigenmodes Y n (r) with 
energies below some classical-field cutoff e cut (the choice 
of which is discussed in Appendix [Ah, and a complemen- 
tary incoherent region I = {n : e cut < e n < E max }. 
Introducing the projector 



P c {/(r)}=£y n (r)/dr' 
nec J 



F n *(r')/(r'), (4) 



onto the coherent region C, we define a C-region field 
operator 



7/>c(r) = 7>c{* L (r)} = 



nec 



,Y n (v). 



(5) 



In the SPGPE formalism, the C-region field operator 
V>c( r ) is treated in an open-systems approach, and the 
complementary I region of the field is regarded as a ther- 
mal and diffusive bath to which the C region is coupled. 
The resulting master equation for the C-region density 
operator corresponding to ipc (r) is simplified by a high- 
temperature approximation and (after neglecting terms 
which do not affect the equilibrium properties of the 
system [44]) is mapped, using standard techniques [49] . 
onto a stochastic field equation in the Wigner represen- 
tation J5UJ EI] for a classical field 



V>cM)= ^a n (t)Y n (r). 



nec 



The resulting equation of motion 



# c (r,i) =V C { - ^C c i> c {r,t)dt 



(6) 



(7) 



+ -^[ f i-£c]ipc(v,t)dt + dW 7 (r,t) j>, 



is termed the simple growth SPGPE [321144] . 

The growth rate 7 quantifies the strength of thermal 
and diffusive damping of the C-region field ^c(r,t) by 
the high-energy bath of atoms in I, and dW 1 {r,t) is a 
complex stochastic noise term associated with this damp- 
ing, which satisfies 



{dW:{T,t)dWJv',t)) = 2 7 8 c (r,r')dt, 



(8) 



where <5c(r,r') = ^2 n(£C Y n (r)Y* (r') acts as a Dirac 
delta function within the C region. The Hamiltonian 



evolution operator for the C region, Cc , is defined by its 
action on the C-region field: 

£ c yj c (r,t) = (k P + WcM)| 2 ^c(r,i)- (9) 

Neglecting all but the first term on the right-hand 
side (RHS) of Eq. (TtJ) we obtain the projected Gross- 
Pitaevskii equation (PGPE) [52"H54"] . The second term 
on the RHS of Eq. ^ is dissipative and, in general, in- 
duces changes in the population Nc — Jdr |-0c(r)| 2 and 
energy E c = Jrfr^(r)[ff sp + (C/ /2)|^ c (r)| 2 ]Vc(r) of 
the C-region field. Within this term, Cc can be thought 
of as extracting the effective (local) chemical potential of 
the classical field; neglecting for simplicity the phase of 
the field, the local field amplitude ipc( r ) therefore grows 
where £c^c(r) is smaller than /i-0c(r), and vice versa. 
The complex noise term dW n (r, t) reflects the stochastic 
nature of the dissipation, which results physically from 
the random scattering of atoms into and out of the C 
region. 

The growth rate 7 can be derived systematically from 
the thermodynamic parameters of the bath and the 
choice of energy cutoff e cut [H]. However, in the present 
study the precise value of 7 is unimportant, as we are only 
concerned with the equilibrium properties of the system, 
and not the detailed nonequilibrium dynamics of its re- 
laxation. We thus choose a value for 7 on the basis of 
numerical expediency (see Sec. Ill A). 



B. Calculation of observables 

The noise and damping terms in the SPGPE [Eq. ([7|] 
serve to drive trajectories of the classical field to a grand- 
canonical equilibrium distribution consistent with the im- 
posed (thermal bath) temperature T and chemical po- 
tential fi [32 . The field undergoes a period of non- 
equilibrium dynamical evolution as it thermalizes toward 
equilibrium with the (above-cutoff) I region. The growth 
of the (quasi-)condensate from an evaporatively cooled 
thermal cloud would, in principle, be modeled by start- 
ing with an initial state '0c(r,' : = 0) corresponding to 
a high-temperature, noncondensed field [Ml [SSJ. How- 
ever, as our interest here is in the equilibrium config- 
urations of the system, we choose for our initial state 
the ground state of the Gross-Pitaevskii equation, which 
we obtain for each considered chemical potential /1 by 
imaginary-time evolution. In this way we avoid the spon- 
taneous formation of long-lived phase defects during the 
passage of the system from the noncondensed to the 
(quasi-)condensed phase [55H57] . which can significantly 
delay complete thermalization [55] . 

Once equilibrium is established, we characterize the 
equilibrium state of the system by calculating correlation 
functions of the classical field; i.e., averages of functionals 
.F[V>c(r)] of the field over the equilibrium distribution of 
field configurations [54 . The I region itself is modeled 
in a semiclassical Hartree-Fock approximation [5S1 [5j|] 
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(see Appendix [b]), from which we can infer the total field 
density n(r) = Uc(r) + n.i(r), and the total atom num- 
ber N — J drn(r) = Ac + N\. The classical correlation 
functions of V'c(r) are the classical-field analogs of quan- 
tum correlation functions of the Bose field ipc(v) [54] - 
and we interpret them as estimates of the corresponding 
quantum correlation functions (i.e., we neglect the for- 
mal commutator corrections of the Wigner theory |32j . 
which is equivalent to neglecting quantum fluctuations). 
In practice, we substitute time averages of a single tra- 
jectory ipc( r ,t) for averages over the grand-canonical en- 
semble: 



Wc(r)])4E*(M,)], (10) 

S 3=1 

following the ergodic interpretation of the formally mi- 
crocanonical classical-field methods [50] [5Tj . 



1. Coherence and (quasi-) condensation 

In this article we consider only equal-time correlations 
of the classical field ipe- The expectation values of all 
one-body observables in the field at equilibrium are en- 
coded by the first-order coherence function [BS] [S3] 



G«(r 



(^(r)Vc(r')). 



(11) 



The local (r' = r) first-order coherence function yields 
the mean density of atoms in the C region, ric(r) = 
G^\r,r), while the off-diagonal (r' =^ r) elements 
G^\r,r') depend additionally on the coherence of the 
phase <j>(r) of the classical field [defined by V'c(r) = 
|V'c( r )|e I '^ r ' ) ] between positions r and r'. The matrix 
G*- 1 ) is Hermitian, and can therefore be diagonalized to 
obtain a complete basis of eigenvectors with real eigen- 
values. Transposing the Penrose-Onsager definition [16] 
of BEC to the classical-field description, we identify the 
largest of these eigenvalues as the condensate popula- 
tion No, and the corresponding eigenvector ipo(r) as the 
(unit-normalized) condensate orbital. This identification 
is supported ex post facto by a consideration of more 



general field correlations, as we discuss in Sec. Ill 



The inhomogeneous system we investigate may exhibit 
more general quasicondensate behavior, associated with 
slow decays of phase coherence across the sample (26] 
|2"51 [50] , We will therefore make use of the normalized 
first-order coherence function 



5 W(r,r') 



GW(r,r') 
yJn c {v)n c (Y') ' 



(12) 



to characterize the spatial decay of phase coherence in 
the system. 

Another important characterization of the field fluc- 
tuations is given by the normalized local second-order 



coherence function 

S (2) (r,r) 



(iV'c(r)l 4 ) 

(hMr)| 2 ) 2: 



(13) 



which is directly related to the density variance: 
Var{nc(r)} = (g( 2 )(r,r) — l)nc 2 (r) [64], and there- 
fore provides a measure of density fluctuations. In 
the limiting case of a purely thermal (chaotic) field, 
</ 2 )(r,r) = 2, whereas g^(r,r) = 1 for a perfectly co- 
herent field [541 165] . 

In a homogeneous system, a quasicondensate is a com- 
ponent of the field that undergoes large point-to-point 
phase fluctuations, but in which density fluctuations are 
suppressed [H] [15] , The quasicondensate density can 
therefore be estimated by considering the extent to which 
fluctuations of the field fail to be Gaussian [55] [ST]- Keep- 
ing in mind that such identifications bear in general no 
a priori relation to the structure of nonlocal phase cor- 
relations in the field, they nevertheless provide a useful 
characterization of the quasicondensate [66, 67] ■ Gener- 
alizing to the inhomogeneous case, the quantity [63j 



(14) 



n QC (r) = J2n c 2 (r)-G( 2 )(r,r), 



provides a useful measure of the field statistics in exper- 
imentally relevant systems 36, 68J, and is often simply 
referred to as the quasicondensate density — a terminol- 
ogy we will also adopt in this article. Following from this 
definition we define the total quasicondensate population 
^QC = /*tiQc(r). 



2. Anomalous correlations 

We further characterize the physics of our system by 
analyzing the so-called anomalous correlations of the field 
fluctuations. In contrast to traditional mean-field theo- 
ries of BEC (see, e.g., Ref. [12]), our approach does not 
assume a priori a particular fixed value for the conden- 
sate phase, and hence preserves U(l) phase symmetry, 
in the sense that the grand-canonical ensemble sampled 
by the classical field is invariant under global phase ro- 
tations. 

Following Ref. [69] (see also Refs. [37] [JDHJJ]), we de- 
fine the fluctuation field 



AM) 



a* (t) 



^a*(t)a (t) 



#c(r,*), 



(15) 



where ao(t) = Jdr (y9g(r)-0c(r, t) is the classical-field am- 
plitude of the (unit-normalized) condensate mode ipo [73j , 
and 



Sip c (r,t) = i/ic(r,t) - a (t)tp (r), 



(16) 



is the component of the classical field ipc(^> t) orthogonal 
to the condensate. Introducing A(r) allows us to calcu- 
late anomalous moments of the non-condensed compo- 
nent of the field, and we will consider in particular the 
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so-called anomalous thermal density (or anomalous av- 
erage) n(r) = (A(r)A(r)) [63]. Whereas the quantity 
(dtpc(' r )S'>Pc( r )) vanishes in the grand-canonical ensem- 
ble sampled by the SPGPE, which is symmetric with 
respect to global phase rotations, k(t) can acquire a 
nonzero value in this ensemble, giving a measure of "pair- 
ing" correlations induced in the noncondensed compo- 
nent of the field by the condensate [H] . Physically, these 
correlations arise due to the coherent scattering of pairs 
of atoms out of the condensate [46] , and the correspond- 
ing time-reversed processes. At equilibrium the rates of 
forward and reverse scattering must balance, and so we 
expect «(r) to be purely real, relative to a real condensate 
orbital ESim 



III. RESULTS 
A. Physical parameters 

We consider a finite-temperature cloud of 23 Na atoms 
in an elongated (cigar-shaped) harmonic potential with 
longitudinal and transverse trapping frequencies uj z = 
2-7T x 5 Hz and u x — u y = uj± = 2n x 250 Hz, respectively. 
Results for some physical quantities will be specified in 
terms of the (long-axis) oscillator length z = ^h/mu) z 
and oscillator energy tuo z . We focus primarily on the 
dependence of the field correlations on the system tem- 
perature at constant total (C-region plus I-region) atom 
number N = 2 x 10 5 . 

As the system we study is fundamentally three- 
dimensional, we expect the transition temperature of 
the corresponding ideal Bose gas to provide a useful 
point of comparison (in contrast to inherently one- and 
two-dimensional systems, in which interactions can eas- 
ily render the corresponding ideal-gas models irrele- 
vant [7U [75]). We calculate the critical temperature 
for the Bose gas in the elongated potential in the non- 
interacting limit 



ftuj 
k B 



— V 

cm) 



1 C(2) Urn 

2 C(3) Q 



(17) 



where C is the Riemann zeta function, ui = (cJxUJybJz) 1 ' 3 



and id- 



w z )/3. The second term on 



in Appendix [X] Characterizations of the system at con- 
stant chemical potential /i, and at constant temperature 
T, are reported in Appendix [C] 

As our classical-field method is grand canonical, in or- 
der to perform calculations at fixed total atom number N, 
we must first calculate the appropriate chemical potential 
fi(N,T). This is achieved using a semiclassical Hartree- 
Fock description of the above-cutoff atoms ( Appendix[B| . 
For simulations performed at fixed \i or fixed T (Ap- 
pendix [C]), we simply calculate the above-cutoff density 
and total atom number after performing our classical- 
field simulations. We choose values of the growth rate 
7 corresponding to a fixed dimensionless growth coef- 
ficient -i/k B T = 40ma 2 k B T/irh 3 in all our simulations 
(cf. Rcfs. [32, 44 J. We find for our choice of parameters 
that thermalization in each simulation occurs over a time 
scale of approximately 100 ms of system time, though to 
ensure sufficient equilibration we form our time averages 
from 10 equally spaced samples taken over a period of 
100 s, beginning after the first 1 s of evolution. 



B. Condensate and quasicondensate 

In Fig. Hi we show the condensate density tiq(z) — 
N \ip (zz)\^ (dotted blue line) and the quasicondensate 
density uqc(z) = nQc(zz) (dashed blue line), on the 
z axis, at a set of representative temperatures. We also 
plot the anomalous thermal density k(z) — k(zz) (dot- 
dashed red line) , and the total density of the Bose field 
n(z) = nc(zz) + m(zz) (solid red line) [see Eq. ( |B1[ )]. 
At the lowest temperatures considered [Figs. [l|a)-[ljc)] 
we observe significant condensate density, and associated 
anomalous thermal density. We find that the anomalous 
density is nonzero only in the central region of the trap 
where the condensate exists, and its magnitude exhibits 
a pronounced dip at the center of the trap, consistent 
with the results of previous works [37 ] l69 [ ITO ] I78H8T] . 

We note that even in this very low temperature regime, 
the quasicondensate density uqc(z) is somewhat larger 
than the condensate density uq(z) [S2]. However, this 
is also a feature of condensation in a less anisotropic 
(oblate) 3D trap, as observed in Ref. [53]. In particu- 
lar, the small excess proportion of quasicondensate here 
is largely attributable to the presence of the anomalous 
average; i.e., it is mostly accounted for by the (nega- 



correction to the critical temperature due to the trap 



the RHS of Eq. p7| is the (leading-order) finite-size tive) contribution of (v?q) 2 (AA) + c.c. (not shown) to the 

local coherence function G^(r,r) (see Ref. [69 ). In 
the limit that the system exhibits well-defined Bogoli- 
ubov quasiparticles, the anomalous average is a measure 
of the phase-fluctuation-like nature [53] of the lowest- 
lying excitations. The phase-fluctuation character of 
these excitations implies that their thermal population 
contributes little to the density fluctuations of the sam- 
ple, as compared (for example) to dressed single-particle 
(Hartree-Fock) states [84] . This effect can be exhibited in 
standard mean-field calculations that linearize the field 
fluctuations about a well-defined condensate mode 1851. 



anisotropy [76l [77] . For our system Eq. ( 17) gives a tran- 
sition temperature T° ss 174 nK. At temperatures be- 
low T°, the condensate fraction of the ideal gas varies as 
JV /iV« 1- (T/T c ) 3 [77]. 

We simulate the field at a range of temperatures T = 
85 - 180 nK. This range extends from the lowest tem- 
perature at which the validity criteria of the SPGPE can 
be satisfied to just above the ideal-gas estimate of the 
BEC transition temperature [Eq. (17)]. We discuss the 



rationale behind our choice of system parameters further 
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(d) T = 120 nK 



( g ) T = 160 nK 
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FIG. 1: (Color online) Variation of system densities on the 
long (z) axis with temperature: condensate density no(z) 
(dotted blue line), quasicondensate density uqc(z) (dashed 
blue line), anomalous density k(z) (dot-dashed red line), and 
total (C-region plus semiclassical I-region) density n(z) (solid 
red line). 



and is implicit in the BdG density-phase approach of 
Refs. [2j)] [25]. In the sense that the field correlations 
at these low temperatures can be understood in terms 
of a distinct condensate with well-defined quasiparticles, 
they are consistent with true condensation. 

As the temperature is incr eased [Fig. [jjd-f)], the size 
of the quasicondensate progressively decreases, and the 
atomic density is increasingly redistributed to purely 
thermal population in the wings of the cloud. At the 
same time, the condensate density recedes from the qua- 
sicondensate density and the anomalous density corre- 
spondingly becomes smaller. At reasonably high tem- 
peratures [Fig. Qg-h)], meaningful condensation (as ev- 
idenced by the persistent presence of an associated 
anomalous average) remains, though both the conden- 
sate and anomalous average are significantly smaller than 
the quasicondensate. In this regime, the large excess of 
quasicondensate density is no longer accounted for by 
the presence of the anomalous average. In particular, 
the presence of significant quasicondensation outside the 
spatial extent of the condensate indicates that the qua- 
sicondensate is largely composed of higher-order corre- 
lations in the field 69 , reflecting the breakdown of a 
description in terms of Bogoliubov quasiparticle excita- 
tions about a well-defined condensate [IBJ [SB]. In this 
regime the alternative approach of Bogoliubov theory 



in phase and density fluctuations [30] is required. Our 
SPGPE approach, which naturally accounts for all dy- 
namical processes in the field (within the classical-field 
approximation) , is equally applicable to the quasiconden- 
sate regime, the low-temperature true condensate regime, 
and the crossover between the two. At the highest tem- 
perature shown [Fig. [lji)], a small quasicondensate re- 
mains, while the condensate (and the associated anoma- 
lous average) is vanishingly small. 

We note that the progression of field correlations with 
increasing temperature shown here is qualitatively sim- 
ilar to that obtained in classical-field calculations for a 
more standard 3D geometry [BJ|]. In both cases, as the 
temperature of the field is increased, the condensate be- 
comes an increasingly small proportion of the total quasi- 
condensate, and the quasicondensate is increasingly com- 
posed of correlations beyond a simple Gaussian (Hartree- 
Fock-Bogoliubov [JT]) ansatz for field fluctuations about 
a well-defined condensate. The difference between the 
present scenario and the more typical 3D geometry is 
therefore largely a quantitative one: in the present case, 
the divergence of the condensate from the quasiconden- 
sate occurs at significantly lower temperatures, and the 
quasicondensate-dominated regime thus extends over a 
large temperature range [57] . 

In Fig. [2ja) we plot the condensate and quasiconden- 
sate populations (Nq and Nqc, respectively) as frac- 
tions of the total atom number N. The occupancies 
of the condensate and quasicondensate exhibit similar 
trends, smoothly increasing from near zero at high tem- 
peratures toward N as the temperature decreases toward 
zero. Interestingly, the functional dependence of Nq on 
temperature appears to correspond more closely to the 
No/N ~ 1~T/T C scaling of the condensate fraction of an 
ideal gas in one-dimensional harmonic confinement [88] . 
than that of the full 3D geometry (Sec. Ill A I. 

The marked difference between the occupancies of the 
condensate and quasicondensate is illustrated by their ra- 
tio Nq/Nqq, plotted in Fig. f2Vb) (asterisks). This ratio 
is small at temperatures T J> 160 nK, and the system 
is therefore dominated by quasicondensate behavior in 
this regime. At lower temperatures this ratio steadily in- 
creases with decreasing temperature, indicating a gradual 
crossover from a quasicondensate to a true condensate. 



C. Identification of the BEC critical point 

1. Binder cumulant analysis 

In Fig. [5Jb) we also plot the generalized Binder cu- 
mulant Cb = (l a o| 4 )/(| Q! o| 2 ) 2 (squar es), w here ao(t) is 
the condensate amplitude (see Sec. II B). Bczctt and 
Blakie [5j5] introduced this quantity as the natural gener- 
alization of the Binder cumulant [90] of a homogeneous 
system to the harmonically trapped case. In homoge- 
neous models within the 3D XY universality class, the 



Binder cumulant acquires the universal value (Cb ) cr it 
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FIG. 2: (Color online) (a) Dependence of condensate (pluses) 
and quasicondensate (crosses) fractions on system tempera- 
ture. Also shown are the ideal-gas BEC transition tempera- 
ture T° (vertical dotted line) and condensate fraction (dot- 
dashed line), (b) Ratio of condensate and quasicondensate 
populations No/Nqc (asterisks), and generalized Binder cu- 
mulant Cb (squares). The horizontal dashed line indicates 
the critical value (C_B) C rit, and the vertical dashed line indi- 
cates the corresponding critical temperature. The solid line 
is a smooth line of best fit used to identify the intersection 
with (Cb )crit- (c) (Negative of the) integrated anomalous 
density (circles), and ratio of expectation values of the ef- 
fective many-body T matrix and two-body T matrix in the 
condensate (diamonds). 



1.2430 at the critical point associated with the transi- 
tion to long-range order [§T]. This fact was used to 
identify the critical temperature in classical-field sim- 
ulations of the 3D homogeneous Bose gas [92] • More 
recently, we found [35] that an estimate of the critical 
point of a harmonically trapped Bose gas based on the 
condition Cb = {Cb) ait was consistent with an indepen- 
dent estimation based on the suppression of condensate- 
condensate interactions due to many-body effects |93j . 
Here we take the condition Cb = {CB)ait as an estimate 
of the location of the Bose-condensation transition in the 
elongated 3D system, and thereby identify the critical 
temperature T c — 165 nK. Our analysis of the Binder 
cumulant therefore indicates a transition to Bose con- 
densation at a temperature slightly below the appropri- 
ate ideal-gas BEC temperature, consistent with previous 
results for a less anisotropic harmonically trapped 3D 
system [2U1IM]. 



2. Many-body T-matrix analysis 

In Fig. [2jc) we plot (the negative of) the integrated 
anomalous density fdrn(r) (circles), which reaches its 
maximum absolute value at intermediate temperatures, 
consistent with previous studies of finite-temperature 
condensates [37l [69] EE ESHEI]. As we have noted 
(Sec. II B), the anomalous density arises due to so-called 
Bogoliubov processes, in which pairs of atoms scatter 
each other out of the condensate and into the noncon- 
densed modes of the field (and vice versa). These pro- 
cesses yield a correction to the strength of condensate- 
condensate interactions, due to Bose-stimulated scatter- 
ing of colliding condensate atoms through occupied in- 
termediate modes. In mean-field theories, this correction 
is encoded in the replacement of the s-wave (contact- 
potential [SS]) two-body T matrix T 2 b(r, r') = UoS(r— r') 
by an approximate many-body T matrix [3S1 l?8Tl8T)l l9"Bl - 

m 



r mb (r,r') = [/o 1 



«(r) 



<J(r- 



(18) 



To quantify the total correction to condensate- 
condensate scattering due to this many-body effect, we 
calculate the matrix element of T m b in the condensate 
orbital 



(T m b)o = ( < /Wo|T' m b|</Wo) 



= U a / dr 1 



K ( T ) 
NofUr) 



Mr)| z 



(19) 



and compare it with the corresponding matrix element 
of the two-body T matrix, (T 2b )o = U fdr |^o0)| 4 . The 
ratio (T m b)o/(72b)o of the two matrix elements is plotted 
in Fig. [2^c) (black diamonds) . 

As is well known, the effect of the many-body processes 
encoded in T m b is to suppress the strength of scattering 
between condensate atoms: In a homogeneous 3D sys- 
tem the effective interaction strength between conden- 
sate atoms vanishes at the BEC phase transition [9"9T - 
1101] . More generally, one finds in inhomogeneous geome- 
tries that the ratio (T m b)o/(l2b)o exhibits a minimum at 
the critical temperature for Bose condensation [69] ITST 
I501I102) . From Fig.[2jc) we observe that the ratio of effec- 
tive T-matrix elements indeed exhibits a minimum value 
close to the critical temperature T c estimated from the 
analysis of the Binder cumulant (Sec. Ill C I. This behav- 
ior strongly suggests that despite the phase-fluctuating 
(quasicondensate) nature of the system at temperatures 
T < T c , the system does indeed exhibit a second-order 
phase transition to condensation, and moreover provides 
an independent validation of the Binder cumulant condi- 
tion Cb = (Cs)crit used to estimate the critical point. 



D. Phase coherence 

To characterize the temperature dependence of the 
equilibrium behavior of the system in more detail, we 



analyze the normalized first-order correlation function 
on the z (long) axis, g^'(z,z') = g^\zz, z'z) : which 
reveals the spatial extent of phase coherence in the sys- 
tem. In Fig. [3] we plot the centrally evaluated coherence 
function g^(0 7 z) (dot-dashed green line) and the sym- 
metrically evaluated coherence function (—z, z) (solid 
green line), along with the condensate density n (z) 
(dotted blue line) and quasicondensate density uqc{z) 
(dashed blue line), for the representative temperatures 
considered in Fig. [I] The condensate mode shape is of 
course determined by the form of g^(r,r') through its 
definition in terms of the Penrose-Onsager criterion. By 
contrast, the density profile uqc(z) of the quasiconden- 
sate, which simply corresponds to the suppression of local 
density fluctuations, is a priori unrelated to the first- 
order coherence of the system (see Sec. II B|. 

Above the critical point [Fig. [3Vi)] , the spatial extent 
of both </ 1; (0, z) and g^ (— z, z) is very narrow, with 
both functions decaying on the length scale of the ther- 
mal de Broglie wavelength At = Hw2tt /mksT ', which in 
this regime is far smaller than the axial harmonic oscil- 
lator length zq = yjhjmujz 62J. Correspondingly, the 
condensate itself is very small, whereas the quasiconden- 
sate is still finite across a comparatively large spatial ex- 
tent, indicating a regime of density-fluctuation suppres- 
sion without any significant one-body coherence. 

At and slightly below the critical temperature 
[Figs. [3jli) and[3^g), respectively] we observe markedly 
broader profiles for both g^(0, z) and g^(— z, z). More- 
over, the symmetrically evaluated correlation g^'(—z,z) 
exhibits a roughly exponential decay with z, as expected 
in the quasicondensate regime [23 1103] . The profiles of 
both these correlation functions remain far narrower than 
tiqc (z) , indicating that the system is well within the qua- 
sicondensate regime, consistent with the small magnitude 
of the condensate density 710(2:) relative to that of the 
quasicondensate tiqc (z) . 

At lower temperatures [Figs.[3ja)-[3]jf)] we see evidence 
of a crossover from a quasicondensate to a true con- 
densate: The profiles of g^(0,z) and g^(— z, z) both 
broaden and become increasingly comparable in width to 
the quasicondensate density riQc (z) with decreasing tem- 
perature. Furthermore, the shape of g^ (— z, z) departs 
from exponential decay and comes to more closely resem- 
ble Gaussian decay at these lower temperatures |103j . 
particularly in Fig. [3^a-d), and the condensate density 
no(z) approaches Uqc(z) with decreasing temperature. 



E. Identification of the 
condensate/quasicondensate crossover temperature 

To quantitatively characterize the crossover from qua- 
sicondensate to true condensate, we fit the function 

f(z) = e -[( 1 -«)( 2z / i *)+«( 2z / i *) 2 ] (20) 

to the correlation function g^\— z, z), where is the 
phase-coherence length, and £ € [0, 1] is a second fitting 
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FIG. 3: (Color online) First-order correlation functions 
g^(— z, z) (solid green) and g^\0, z) (dot-dashed green), and 
densities of the condensate 72,0(2) (dotted blue) and quasicon- 
densate nQc(z) (dashed blue). 



parameter that controls the functional form of f(z). The 
function f(z) exhibits purely exponential decay in the 
limit £ = 0, characteristic of g^'{— z, z) in the quasicon- 
densate regime; and purely Gaussian decay in the limit 
£ = 1, characteristic of g^ (— z, z) in the true condensate 
regime [5(3 1103] . Using this fitting function has the ad- 
vantage of taking into consideration both the width and 
shape of g^~> (-Z, z) in providing an estimate of the coher- 
ence length L<j, [35]. Moreover, it allows us to identify a 
qualitative change in the correlation function [3UJ H03] as 
a signature of the crossover, and in this article we define 
the characteristic temperature of the crossover as that at 
which £ = 0.5 (cf. Ref. [I)]). 

In Fig. |4|a) we plot the temperature dependence of £ 
(empty stars), and observe that the characteristic value 
£ = 0.5 of the crossover (horizontal dashed line) is 
reached at a temperature = 135 nK (vertical dot- 
dashed line). We observe a gradual change in £ through 
the crossover, from £ = above T c (vertical dashed line) 
to £ = 1 far below T c , consistent with the gradual change 
in the shape of g^(— z, z) shown in Fig. pi 

In Fig. [4](b) we plot the temperature dependence of 
(empty circles), along with the estimated Thomas- Fermi 
lengths of the condensate (pluses) and quasiconden- 
sate L^ C ^ (crosses). The two Thomas-Fermi lengths are 
obtained by performing fits of the condensate and qua- 
sicondensate column densities, 7To(y, z) = fdxno(r) and 
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FIG. 4: (Color online) (a) Fitting parameter £, which charac- 
terizes the functional form of g^ 1 ' (— z, z) (see text), in SPGPE 
(empty stars) and BdG |28| (filled stars) calculations. The 
horizontal dashed line indicates the condition £ = 0.5, by 
which we define the crossover temperature. The solid lines are 
smooth lines of best fit used to identify the intersection with 
£ = 0.5. (b) Phase coherence length obtained from SPGPE 
(empty circles) and BdG [28 s (filled circles) calculations, and 
fitted Thomas-Fermi lengths Li ' (pluses) and Z/!^ ' (crosses) 
of the condensate and quasicondensate, respectively. Vertical 
lines through both panels indicate the crossover temperature 
T$ obtained from SPGPE (vertical dot-dashed line) and BdG 
[28| (vertical dotted line) calculations, and the critical point 
identified using Cb (vertical dashed line). 



nQc(y, z) = JdxnQc(r), respectively, to the Thomas- 
Fermi column density 



Utf(2/, z) = 7t,tf(0, 0)max|o, ^1 



y 

Li 



z 

Li 



'}, (21) 



obtained by integrating Eq. (B3) over x. Here we regard 
the peak column density titf(0, 0) , transverse Thomas- 
Fermi length Lj_, and axial Thomas- Fermi length L z as 
independent fitting parameters [104| . Previous studies 
of this system loosely distinguish the quasicondensate 
and true condensate regimes by the size of the coher- 
ence length as compared to the Thomas-Fermi length of 
the quasicondensate [251 [35]. Indeed we find here that 
at temperatures in the range < T < T c , the coher- 
ence length Z/0 is smaller than the axial length of the 

quasicondensate L^ C \ whereas the coherence length 

increasingly exceeds as the temperature decreases 

below T,/,. We note in particular that equality of the two 

lengths (L^ — L^ C ^) appears to coincide almost exactly 
with the crossover temperature as defined by the condi- 
tion f = 0.5. 

We now compare the crossover temperature obtained 
from our SPGPE simulations with the predictions of the 
BdG theory of Petrov et al. EE]. Th e authors of Ref. [2S] 



derive an approximate expression 
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for the variance of the phase difference between ax- 
ial points z and z', where n(No) and L z (Nq) are the 
Thomas-Fermi chemical potential and associated axial 
Thomas-Fermi length, respectively, of a pure condensate 
of population No [we use the ideal-gas estimate of the 
condensate population N (T) — see Sec. |III A] , and the 
Pj are Jacobi polynomials. Using Eq. (pij) we calcu- 
late the symmetrically evaluated coherence function 

5 «(-z,z) = e-^-^>] 2 >, (23) 

for each temperature T we simulated with the SPGPE. 
We then proceed to determine values of £ and by 
fitting g( x '{— z,z) to Eq. (20), and plot these quantities 
in Fig. [4] (filled stars and filled circles, respectively) where 
they may be compared to our SPGPE results. 

We find that the BdG model predicts the crossover cri- 
terion £ = 0.5 to occur at a temperature = 144 nK 
(vertical dotted line) , which we note lies somewhat above 
the temperature at which it occurs in our SPGPE calcu- 
lations (vertical dot-dashed line) . Proukakis H05| has 
shown that in true one-dimensional systems, the density 
fluctuations neglected in the treatment of Refs. [26l [28] 
can erode phase coherence in the system, and thus push 
the crossover to lower temperatures. We note that in ad- 
dition, our 3D classical-field model includes fluctuations 
of the field along the transverse (tight) axes of the trap, 
which are neglected in the model of Ref. [28:, and may 
further reduce the phase coherence in the system. Nev- 
ertheless, our quantitative results show that the analy- 
sis of Petrov et al. [28 j captures the essential physics of 
the quasi-condensate regime of degenerate Bose gases in 
elongated 3D traps. 



IV. CONCLUSIONS 

We have investigated the correlations of a weakly in- 
teracting Bose gas in an elongated 3D harmonic trap- 
ping geometry, using a grand-canonical classical-field 
method. Our investigations spanned temperatures rang- 
ing from the low-temperature true-condensate regime, 
through the crossover to the phase-fluctuating quasicon- 
densate regime, to the normal phase of the gas at high 
temperatures. We characterized the onset of condensa- 
tion in the system using two independent measures: the 
Binder cumulant quantifying number fluctuations in the 
Penrose-Onsager condensate orbital, and the suppression 
of the effective two-body interaction strength by many- 
body processes, as encoded by the anomalous thermal 



10 



density of the field. We found that both measures in- 
dicated that the transition to condensation occurs at a 
temperature slightly below the ideal-gas critical temper- 
ature for this geometry, as is the case in more standard 
3D geometries [2"9" ] l9l ] . 

However, our results show that the system remains in a 
strongly phase-fluctuating quasicondensate regime until 
significantly lower temperatures, as previously predicted 
for such an anisotropic geometry [28] . We explained that 
the quasicondensate phase should fundamentally be un- 
derstood in terms of partial Bose condensation of the 
field, but that the large phase fluctuations in this regime 
imply that the noncondensed component of the gas can- 
not be understood in terms of Gaussian (Hartree-Fock- 
Bogoliubov) fluctuations. 

We identified a temperature characteristic of the 
condensate-quasicondensate crossover, based on a quali- 
tative change in the functional form of the non-local first- 
order coherence function, and found that this qualitative 
change occurs at roughly the same temperature at which 
the phase-coherence length equals the axial length of the 
quasicondensate. The crossover temperature we find is 
somewhat lower than that predicted by an approximate 
Bogoliubov-de Gennes model of phase fluctuations in the 
system [35] , consistent with the expectation that density 
fluctuations of the field have an additional deleterious 
effect on phase coherence in the system [351 H05] . 

The pseudo-low-dimensional system we have consid- 
ered is conceptually simpler than true low-dimensional 
systems, in that the appearance of the condensate, which 
here underlies the quasicondensate regime of the gas, 
is fundamentally due to a proper thermodynamic phase 
transition; i.e., it is expected to persist in the thermo- 
dynamic limit. By contrast, in true low-dimensional 
systems, a significantly occupied orbital that can mean- 
ingfully be identified as a condensate may simply occur 
due to finite-size effects. Nevertheless, the tools used in 
our characterization of the role of the condensate in the 
phase-fluctuating quasicondensate regime of this compar- 
atively straightforward geometry offer to help elucidate 
the role of the Penrose- Onsager "condensate" in truly 
(quasi-)low-dimensional inhomogeneous Bose systems. 
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Appendix A: Simulation parameters 

In choosing simulation parameters we must first en- 
sure that our system satisfies the physical requirements 
to exhibit elongated-3D quasicondensate behavior over 
an appreciable temperature regime. Using the estimates 
of Petrov et al. [28] for T c and T in the limit N « N, 
we obtain the approximate scaling 



fH 3 N 4 



\m 3 a 6 uj^ 



(Al) 



It can be seen that this ratio scales weakly with atom 
number, mass, and scattering length, and that its 
strongest dependence is on the frequencies wj_ and u> z 
of the trapping potential. The (relative) temperature 
range spanned by the quasicondensate regime therefore 
increases with increasing aspect ratio luj_/uj z . However, 
our choice of aspect ratio is limited by the require- 
ments that uj± is not too large, in order that the sys- 
tem remain three-dimensional at the crossover (i.e., that 
fc^T^, fj, > hu>±), and that the value of u> z is experimen- 
tally reasonable. 

We must also take care to choose parameters such that 
the validity conditions of our classical-field methodol- 
ogy GUI can be satisfied. The first such condition is that 
the cutoff energy e cut is high enough that all eigenmodes 
of the single-particle Hamiltonian that are strongly cou- 
pled to one another by interactions (i.e., those which con- 
tribute to the quasicondensate density) are included in 
the C region. We therefore require that the cutoff en- 
ergy satisfies 



(ecut - e )/U n c (0) > 1, 



(A2) 



where eo is the ground-state energy of the single-particle 
Hamiltonian [Eq. and ric(0) is the central (peak) 
density of the C-region field (cf. Ref. [68] )■ We note 
that for a 3D system with fcgT, /j, > huj±, this condition 
automatically implies that the classical field will span 
multiple modes in the transverse dimensions. 

The second c-field validity criterion is that all modes in 
the C region have mean occupations > 1. This condition 
is, in general, at variance with Eq. ( |A2[ ), as raising the en- 
ergy cutoff e cu t to accommodate the effects of interactions 
introduces progressively higher energy — and therefore 
increasingly sparsely populated — modes into the C re- 
gion. Simultaneously satisfying both of these classical- 
field conditions down to low temperatures T < X^, while 
maintaining quasicondensate behavior over a reasonably 
large temperature range is more readily achieved in sys- 
tems with a relatively small dimensionless interaction en- 
ergy Uq = 4iray/ mujz/h, and for this reason we have cho- 
sen to simulate 23 Na atoms. 

For a system of N = 2 x 10 5 atoms of 23 Na, in a trap 
with frequencies (luj_,uj z ) — 27r x (250,5) Hz, we obtain 
in our simulations a temperature ratio of T^jT c ss 0.82 
(see Sees. IIIC and HIE I. These parameters allow us to 



properly satisfy the classical-field validity conditions over 
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FIG. 5: (Color online) Energy cutoff e C ut (circles) and ratio (e cu t — £o)/Uonc(0) of cutoff and interaction energies (triangles) 
for simulations with (a) varying T at constant N, (b) varying T at constant fi, and (c) varying fj, at constant T. 
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FIG. 6: (Color online) C-region population iVc (pluses), and 
semiclassical Hartree-Fock calculations of total (C-region plus 
I-region) atom number N (/i, T) (crosses) and chemical poten- 
tial fj,(N,T) (triangles). Quantities are shown for (a) varying 
T at constant N, (b) varying T at constant fi, and (c) varying 
/i at constant T. 



the temperature range of interest: For each set of param- 
eters we simulated, we chose the value of e cu t (circles in 
Fig.p5| such that rij = (|<x,| 2 ) > 2 for each mode Yj(r) of 
the classical field. The corresponding values of the ratio 
(ecut — £o)/EA) n c(0) > 3/2 (triangles in Fig. [5]) are gener- 
ally sufficient to encompass the quasicondensate density 
within the C region. 
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where the semiclassical energy is 



e(r , p) = f - + V cxt (r) + 2U [n c (r) + m (r)] , (B2) 
2m 

and e(r,p m i n ) = e cut ensures that the integration over 
momentum is restricted to the I region. The semi- 
classical energy includes the mean-field potential from 
both the C-region density nc(r) and the I-region den- 
sity ni(r) [23]. Equation (Bl) must therefore be solved 



iteratively to obtain a self-consistent I-rcgion density pro- 
file for specified values of [i and T. Spatial integration 
of the total density then gives the total atom number, 
N(/i, T) = Jdr [nc(r) + ni(r)]. This can be inverted nu- 
merically to find fJ,(N, T), though doing so requires «c(r) 
to be recalculated via SPGPE simulation and ergodic av- 
eraging after each iteration of fi. In practice, however, 
only a single iteration is required to obtain a value of 
/j.(N,T) that yields a total atom number within 1% of 
the target value N, provided the initial guess for /j, is suf- 
ficiently accurate. We obtain such a guess by numerically 
inverting an approximate form of JV(/k, T), in which we 
replace nc(r) by the Thomas- Fermi condensate density, 

no(r) = max{0, - V ext (r)) /U Q } , (B3) 

and use the condition e(r,p min ) = /i to calculate n t (r). 



Appendix B: Semiclassical Hartree-Fock 

We calculate equilibrium values of N([i,T) and 
/j,(N, T) (plotted in Fig. [6]) using semiclassical Hartree- 
Fock theory as outlined in Refs. [84], IT06HTTT] . For a 
given choice of \i and T, the density of the thermal cloud 
ni(r) is calculated by integrating the Bose-Einstein dis- 



Appendix C: Fixed chemical potential and fixed 
temperature results 

In Fig.[7]we present, in addition to the behavior of the 
system at constant particle number N = 2 x 10 5 over the 
temperature range T — 85 - 180 nK discussed in Sec. [ill 



two additional sets of results for the field correlations: 
those obtained over a range of temperatures T = 50 - 
250 nK at a constant chemical potential /i = 90hu) z , and 
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FIG. 7: (Color online) Properties of the system evaluated at varying T and fixed N (a,d,g,j,m,p), varying T and fixed fj, 
(b,e,h,k,n,q), and varying /i and fixed T (c,f,i,l,o,r): (a-c) Temperature T (squares) and chemical potential /i (triangles), (d-f) 
Condensate fraction Nq/N (pluses) and quasicondensate fraction Nqc/N (crosses). For reference the condensate fraction of 
the ideal Bose gas, and the ideal-gas critical temperature (chemical potential) are also indicated (dot-dashed line and vertical 
dotted line, respectively) [TT2 . (g-i) Ratio JVo /Nqc of condensate and quasicondensate populations (asterisks) and generalized 
Binder cumulant Cb (squares). The horizontal dashed line represents the critical value (Cs) cr it for the Binder cumulant, and 
the vertical dashed line indicates the inferred estimate of the critical temperature. The solid line is a smooth line of best fit used 
to identify the intersection with (Cs) C rit- (j-1) (Negative of the) integrated anomalous thermal density (circles) and effective 
suppression of scattering between condensate atoms due to associated processes (diamonds), (m-o) Fitting parameter £, which 
characterizes the functional form of g^\—z, z) (see text), for SPGPE (empty stars) and BdG [25] (filled stars) calculations. The 
horizontal dashed line indicates the condition £ = 0.5, by which the crossover is characterized for SPGPE (vertical dot-dashed 
line) and BdG [25] (vertical dotted line) calculations. The solid lines are smooth lines of best fit used to identify the intersection 
with £ = 0.5. (p-r) Phase coherence length for SPGPE (empty circles) and BdG [28] (filled circles) calculations, alongside 
fitted Thomas-Fermi lengths Z/i°' (pluses) and L^ ^ (crosses) of the condensate and quasicondensate, respectively. 
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those obtained over a range of chemical potentials [i = 30 
- 230fi.oj z at a constant temperature T = 135 nK. 

Figs.[7]Ja)-[7jc) show the system temperature (squares) 
as a fraction of the ideal gas critical temperature Eq. ( 17) 



[which itself varies with the varying total atom number 
in the data sets corresponding to Fig. Rb) and (c)], and 
the system chemical potential (triangles). Figs.[7](d)-[7]Jf) 
show the condensate fraction (pluses) and quasiconden- 
sate fraction (crosses), together with the ideal-gas result 
(dot-dashed line) for the condensate fraction at the cor- 
responding values of N and T. In all three cases, the con- 
densate and quasicondensate exhibit similar trends, in- 
creasing smoothly with decreasing T (increasing /i) from 
near zero at the highest temperatures (smallest chem- 
ical potentials) considered. Moreover, the quasiconden- 
sate is in general significantly larger than the condensate, 
as shown explicitly by the ratio Nq/Nqc (asterisks) in 
Figs. |7^g)-[7jl). The data also consistently show that Nq 
and Nqq are both suppressed below the corresponding 
ideal-gas predictions for the condensate occupation. 

In Figs. Rg)-[7][i) we plot, as in Fig.^b), the Binder 
cumulant Cb (squares) and observe that the estimate of 
the critical point from the criterion Cb — {Cb) crit um ~ 
formly indicates a lowering of the critical temperature, or 
correspondingly, a raising of the critical chemical poten- 
tial (or density) |104j . This estimate of the critical point 



(vertical dashed line) is seen in all three cases to be rea- 
sonably consistent with the maximal suppression of the 
effective many-body T matrix indicated in Figs.[7]j)-[7]J1) 
(diamonds). Furthermore, we note that the condensate- 
to-quasicondensate population ratio Nq/Nqc [asterisks 
in Figs. [7](g)-[7](i)] assumes a value N /Nqc ~ 0.3 at the 
identifiea critical point in all three data sets. 

In Figs. [7Fm)-[7]jo) we plot the fitting parameter £ 
(empty stars), which c hara cterizes the functional form 
of g^{—z,z) [see Eq. (20)], and identify the crossover 



temperature or chemical potential (vertical dot-dashed 
line) according to the condition £ = 0.5 (horizontal 
dashed line). In all three data sets, our results indicate a 
crossover to true condensation at a temperature (chem- 
ical potential) significantly lower (higher) than both the 
estimated critical point (vertical dashed line) and the 
crossover temperature obtained from the BdG model of 
Petrov et al. [28] (corresponding £ values and crossover 
temperature indicated by filled stars and vertical dotted 
line, respectively). Moreover, we see in Figs. |7£p)-[7[r) 
that the intersection of the phase-coherence length 
(empty circles) and Thomas-Fermi length L^ -* of the 
quasicondensate |104j (crosses) coincides almost exactly 
with the crossover temperature (vertical dot-dashed line) 
in each of the three cases considered. 
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